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ABSTRACT 



Context. Small, bright stellar disks with scale lengths of a few tens of parsec are known to reside in the center of galaxies. They 
are believed to have formed in a dissipational process as the end result of star formation in gas either accreted during a merging (or 
acquisition) event or piled up by the secular evolution of a nuclear bar. Only a few of them have been studied in detail to date. 
Aims. Using archival Hubble Space Telescope (HST) imaging, we investigate the photometric parameters of the nuclear stellar disks 
hosted by three early-type galaxies in the Virgo cluster, NGC 4458, NGC4478, and NGC4570, to constrain the process that forms 
their stars. 

Methods. The central surface brightness, scale length, inclination, and position angle of the nuclear disks were derived by adopting 
the photometric decomposition method introduced by Scorza & Bender and assuming the disks to be infinitesimally thin and expo- 
nential. 

Results. The location, orientation, and size of the nuclear disks is the same in all the images obtained with the Wide Field Planetary 
Camera 2 and Advanced Camera for Surveys and available in the HST Science Archive. The scale length, inclination, and position 
angle of each disk are constant within the errors in the observed U, B, V, and / passbands, independently of their values and the 
properties of the host spheroid. 

Conclusions. We interpret the absence of color gradients in the stellar population of the nuclear disks as the signature that star for- 
mation homogeneously occurred along their length. An inside-out formation scenario is, instead, expected to produce color gradients 
and is therefore ruled out. 

Key words, galaxies: bulges - galaxies: elliptical and lenticular, cD - galaxies: photometry - galaxies: structure 



1. Introduction 

The presence of small and bright stellar disks with scale 
lengths of a few tens of parsec was first discovered 
in the nuclear regions of nearby galaxies about fifteen 
years ago thanks to the subarcsec-res olution capabilities 
of the Hubble Space Telescope (HST, [v an den Bosc h et al 



1994|: iKormendv et al.|[T994l: lEmsellem et a l. 1994; Lauer et al 
' 11995). Afterward, nuclear stellar disks (NSDs) and larger 
scale embedded disks were reported in a large number of 
early-type galaxies as a results of photometric |Lauer et al. | 



X 



1995'; 'Carollo et al.' "1997"; 'Ravindranath et al.' "2001"; "Rest et al.' 
2001; Tranetal. 2001; Ferrarese et al. 2006; Sethetal. 2006; 
Simoes Lopes et al. 2007; Kormendv et al. 200 9j) and spectro- 
scopic surveys (iHallidav et al. 2001; Kuntschner et al.l l2006l: 
iMcDermid et al.ll2006h respectively. Give n that nuclear disk s 
are easiest to detect when nearly edge-on ( Rix & Whitdll99Ql) . 
the observed fraction is consistent with NS Ds be ing a common 
structure in early-type galaxies (iLedo et al .1120101) . At the center 
of spiral galaxies NSDs are also present but these are rel atively 
rare i n late types ff izzella et al. 2002; Falcon-Barros o et al 
20061: iPeletieretal] l2007i: iBalcells et al.l l2007l: iMorelli et al 



2008) 



Photometric parameters of NSDs (i.e., central surface 
brightness, scale-length, axial ratio, and position angle) have 
been derived for nine galaxies spanning a wide range of Hubble 
types, from ellipticals ([Mor elli et al. 2004), to lenticulars 
(Ivan den Bosch et al] 1 19941: iLauer et al.. .1995.: .Kormendv et al.l 



1 19961: Ivan den Bosch et aP 1 19981: IScorza & van den BoschI 
'1998*), and to early-type spirals ( Pizzella et al.ll2002b . They are 
characterized by a smaller scale length (10-30 pc) and higher 
central face-on surface brightness (15-19 mag arcsec"^ in the 
V band) than the embedded stellar disks observed in ellipticals 
and lenticulars as well as to the large kpc - scale disks of lentic- 
ulars and spirals (Ivan den BoschI 119981 : Pizzella et 4] I2QQ2L 
iMorelli e t al. 20041). In particular, the embedded stellar disks 
represent the intermediary between the sequence of elliptical 
galaxies and that of disk galaxies. They have been detected in 
both disky elliptical and lenticular galaxies, have scale lengths 
of 0.1 - 1 kpc, and central face-on surface brightnesses of 
18-21 mag arcsec"^ in the V band. Their angular momentum 
is parallel to that of their host spheroid. This demonstrates that 
the embedded disks are not the result of accretion or merger 
events but are likely to be primordial, in the same sense as the 
main disks of lenticular and spiral galaxies ("Scorza & B ended 
fT995l). The smooth variation in the scale parameters from 
the nuclear to the embedded and main stellar galactic disks 
provides observational support that galaxy properties vary 
with cont inuity along a sequence of increasing disk-to-bulge 
ratio (see 'Kormendy & Bender 1996, and references therein). 
Therefore, unveiling the formation scenario of NSDs may also 
improve our understanding of galaxy formation and evolution. 

In the current picture, NSDs are believed to have formed 
in a dissipational process as the end result of star formation in 
gas either accreted in a merging (or acquisition) event (the most 
clearcut case is NGC 4698. .Bertola et al.lll999l: iPizzella et all 



2 



L. Morelli et aL: Multiband photometric decomposition of nuclear stellar disks 



l2Q02h or pil ed up by the secular evolution of a nuclear bar (e.g., 
NGC 4570, van den Bosch et al. 1998:' Scorza & van den BoschI 
119981: [van den Bosch & Emsellem 1998(). Each of these scenar- 
ios is likely to be correct for some, but not all, objects. In both 
of them, the gas is efficiently directed toward the galaxy cen- 
ter (e.g., [B arnes & Hernquist 1996; Athanassoula et al. 2005; 
Dotti et alEoOTuEliche-Moral et al.i i2009: Hopkins & Ouataert 



2010h . where it first dissipates and settles onto an equilib- 
rium plane and then forms into stars. A striking example of 
an on going p rocess of dissipational formation is provided by 
NGC 4486A (iKormendv et al.l l2005l) . The NSD of this low- 
luminosity elliptical coexists with its progenitor disk of dust 
and gas. The question about the origin of these dynamically- 
cold components is also related to their formation epoch, i.e., 
whether they built up in the early stages of the galaxy assembly 
or whether they formed later. 



To address this iss ue, Ivan den Bosch et al ] p98) 



Ivan den Bosch & EmsellemI (Il998l) . iMorelli et al.1 f2004), 
and iKrainovic & Jaff'd (l2004 studied the age and metallicity 
content of some NSDs. In some cases, the nuclear disk was 
found to be younger than the stellar population of their surround- 
ing sp heroidal component (e.g., NGC 4486A. iKormendy et al.l 
I2QQ5E . in other cases the formation of the disks occurred at the 
same time as the main bod y of the host galaxy (e.g., NGC 4342, 
Ivan den Bosch et al.|[l998h . 

Characterizing how the star formation proceeded in nuclear 
disks can lead to understanding what triggered their forma- 
tion. For instance, if stars formed more or less at the same rate 
throughout the disk from disk material that accumulated very 
rapidly at the center, as one expects to happen during a merging 
event, we should not find strong radial gradients in the proper- 
ties of the stellar populations. In contrast, color gradients are ex- 
pected if the disk experiences an inside-out star formation, l ike 
the disks of spiral galaxies (e.g.. iMunoz-Mateos et al ] l2007h In 
this paper, we use HST archival images of the nuclei of three 
early-type galaxies in the Virgo cluster, NGC 4458, NGC 4478, 
and NGC 4570 to perform a multiband analysis of the photomet- 
ric parameters of their NSDs. Our aim is to understand whether 
the star formation of the nuclear disk was homogeneous and oc- 
curred simultaneously everywhere along its length. 

2. Sample galaxies 

We searched the HST Science Archive for all NSDs known in lit- 
erature that had been imaged in at least three filters with the Wide 
Field Planetary Camera 2 (WFPC2) and Advanced Camera for 
Surveys (ACS) and had known properties of their stellar popu- 
lation. We found multiband images of the nuclear disks hosted 
in three early-type galaxies of the Virgo cluster, NGC 4458, 
NGC 4478, and NGC 4570 (Table[I]). 

NGC 4 458 and N GC 4 478 are classified as EO- 
1 (Ide Vauc ouleurs etaP jl 99 iL hereafter RC3) and E2 
(Sandage & Tamman^ Il98ll hereafter RSA; RC3), respec- 
tively. They are low-luminosity ellipticals with a ~ -18 
(RC3) at a distance of 12.6 Mpc (iTuUvl I1988L Hq = 100 
km s ~^ Mpc~^). Both galaxies host a kinematically-decoupled 



core ( Hallidave t al. 2001; Emsellem et al. 2007). MorelH et al. 
(I2004 detected NSDs inside both on the same radial scale as the 
core and analyzed their structure and stellar content. Their scale 
length and face-on central surface brightness (Table O were 
derived from WFPC2/F814W images and they fit in the relation 
for NSDs. For NGC 4458, these parameters are typical of NSDs, 
while for NGC 4478 they are in-between those of NSDs and the 
disks of disky ellipticals. The mass of the two disks is of few 



10^ Mq. The central disk of NGC 4458 is counter-rotating. It 
has similar properties to the decoupled cores of bright ellipticals 
(see Bertola & Corsini 1999, for a review). According to 
iMorelli et al.l (120041) . it is as old as (^ ^ 15 Gyr) and richer in 
metals ([Z/H] ^ 0.2) than the rest of the galaxy, with a high 
and approximately constant overabundance ([a/Fe] ^ 0.3). The 
NSD of NGC 4478 is co-rotating with respect to the main body 
of its host galaxy, but it is younger (t ^ 6 Gyr), richer in metals 
([Z/H] ^ 0.4) and less overabundant ([a/Fe] ^ 0.2) than the 
outer regions. The nearly solar a-element abundance indicates a 
prolonged star formation history, typical of an und isturbed disk- 
like, gas-rich, and possibly pre-enriched structure. Morelli et al.l 
( 2004) used the models of sin gle stellar popu l ations with a/Fe 
overabundances developed by iThomas et al. 1 (l2QQ3h to derive 
the properties of the nuclear stellar populations. They analyzed 
the absorption line-strength indices measured along the major 
and minor axes of both galaxies. The masses of the two NSDs 
were estimated by a dopting the mass -to-light ratios calculated 
with the models of iMarastonl (1 19981) from the available ages 
and metallicities. The presence of cold disks at the centers of 
NGC 4458 and NGC 4478 argues for a formation scenario based 
on the dissipational collapse of a gas-rich object. An external 
origin of the g as is claimed to ac count for the counter-rotation 
of NGC 4478 (iMorelfi et al.ll2004l) . 

NGC 4570 is a edge-on disk galaxy classified as SO by 
RSA and RC3. It has ~ -19 (RC3) at a distance of 
12.6 Mp c (iTullvl Il988h. Its nuclear disk was first studied in 
detail by Ivan den Bosch et a n (Il998h . It resides within the in- 
ner cutoff* of the main stellar disk, whose light contribution 
to the galaxy surface brightness is negligible within the inner 
2'^ The inner truncation of the main disk is a general feature 
in multi-disk systen i s host ing a nuclear (Ivan den Bosch et al] 
1998"; "Pizzella et al.' '2002' ) and/or an embedde d disk (e.g. 
Scorza & Bender 1995; Sei fert & Scorzd 1996; E msellem et al.l 
1996) . The photometric parameters (Table O were derived by 



analyzing a WFPC2/F555W im age and are typical of NSDs 
(IScorza & van den BoschI IT998I) . The NSD of NGC 4570 has a 
mass of 8 X 10^ M© (at the adopted distance) as derived from dy- 
namical modelling by van den Bosch & Emsellem (1998). It has 
a diff'erent stellar population than the rest of the bulge, consist- 
ing o f relatively old (t ^ 8 Gyr) and metal-rich ([Fe/H] ^ 0.35) 
stars (Ivan den Bosch et aT]|l998h . These results were derived by 
comparing the colors measured in WFPC2 images and the line- 
strength indices measured in Faint Object Spectrograph spectr a 
with the models of single stellar populations by lWortheyl (Il994l) . 

In addiction, NGC 4570 exhibits two nuclear 
edge-on rings. According to dynamical modelling by 
van den Bosch & Emselle based on the observed 

surface brightness and stellar kinematics, the position of the 
rings (at about and 4'' from the center, respectively) is consis- 
tent with the location of the inner Lindblad and ultra-harmonic 
resonances of a rapidly-rotating nuclear bar of about 500 pc. 
The corotation radius is close to the inner cut-off* of the main 
disk at about 7'' from the center. The morphology, dynamics, 
star formation, and chemical enrichment of the nuclear region 
strongly favor a formation scenario of the NSD driven by the 
secular evolution. The nuclear disk appears to be the end result 
of the star formation that occurred in the gas that was funneled 
to the center by the nuclear bar (van den Bosch & EmsellemI 
119981: IKrainovic & Jaff'dl2004 . The growth of a central mass 
concentration produced by the gas inffow altered the shape of 
the orbits sustaining the bar. A s a consequence, the nu clear 
bar weakened and dissolved (see I Athanassoula et al.ll2005l and 
references therein). 



L. Morelli et aL: Multiband photometric decomposition of nuclear stellar disks 
Table 1. Details about the HST archival images of the sample galaxies. 
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Galaxy 


Camera 


Filter 


Obs. Date 


No. Exp. 


Exp. Time 


Prog. Id. 


P.I. 


(1) 


(2) 


(3) 


(4) 


(5) 


[s] 
(6) 


(7) 


(8) 


NGC 4458 


ACS/WFC 


F475W 


28 Dec 2003 


2 


750 


9401 


R Cote 




WFPC2/PC 


F555W 


18 Dec 1998 


3 


1340 


6587 


D. 0. Richstone 




WFr L^Z/r 




iZ fCD Lyyj 


D 


i iZU 


JJ iZ 


o. ivi. raoer 


NGC 4478 


ACS/HRC 


F330W 


28 jul 2005 


4 


1912 


10435 


R. W. O'connell 




ACSAVFC 


F475W 


09 Jul 2003 


2 


750 


9401 


R Cote 




WFPC2/PC 


F555W 


07 Jul 1997 


4 


1600 


6587 


D. 0. Richstone 




WFPC2/PC 


F814W 


07 Jul 1997 


3 


1600 


6587 


D. 0. Richstone 


NGC 4570 


WFPC2/PC 


F336W 


19 Apr 1996 


5 


3100 


6107 


W. Jaffe 




ACS/WFC 


F475W 


13 Jul 2003 


2 


750 


9401 


R Cote 




WFPC2/PC 


F555W 


19 Apr 1996 


2 


400 


6107 


W. Jaffe 




WFPC2/PC 


F814W 


19 Apr 1996 


2 


460 


6107 


W. Jaffe 



NOTES: Col.(l): Galaxy name. Col.(2): Camera name. Col.(3): Filter name. Col.(4): Observation date. Col.(5): Number of exposures. Col.(6): 
Total exposure time. Col. (7): HST proposal number. Col. (8): Name of the principal investigator. 



3. Data reduction 

The WFPC2 and ACS multiband images of NGC 4458, 
NGC 4478, and NGC 4570 were retrieved from the HST Science 
Data Archive. Their list and observing log are given in Table [T] 

The WFPC2 images were obtained by centering the galaxy 
nucleus on the higher resolution Planetary Camera (PC). The 
PC detector is a Loral CCD with 800 x 800 pixels. The pixel 
size is 15 X 15 yum^. The plate scale of 0'.'046 pixel" ^ yields a 
field of view of about 36'' x 36''. To help in identifying and cor- 
recting cosmic-ray events, diff'erent exposures were taken with 
each filter. The ACS images were taken with both the Wide Field 
Channel (WFC) and the High Resolution Channel (HRC). WFC 
consists of two SITe CCDs with 2048 x 4096 pixels each of size 
15 X 15 yum^. The plate scale is 0'.'049 pixel" ^ and the field of 
view of the combined detectors covers an approximately square 
area of about 202" x 202". HRC detector is a SITe CCDs with 
1024 X 1024 pixels each of size 21x21 yum^ with a plate scale 
of a'028 X a'025 pixel"^. Because of the large, but well char- 
acterized, geometric distortion aff'ecting the ACS, the 29" x 26" 
HRC field of view projects onto the plane of the sky as a rhom- 
boid with X and y axes forming a 84? 2 angle. For both WFPC2 
and ACS images, the telescope was guiding in fine lock, giving 
a typical rms tracking error per exposure of 0'.'005. 

The WFPC2 images were calibrated using the CalWFPC re- 
duction pipeline in IRAlQ maintained by the Space Telescope 
Science Institute. Reduction steps include bias subtraction, dark 
current subtraction, and flat-fielding, as described in detail in 
the WFPC2 instrument and data handbooks ( Baggett et al. 2002; 
[McMaster et al.l i2008). Subsequent analysis was performed us- 
ing IRAF standard tasks. The bad pixels were corrected by 
means of a linear one-dimensional interpolation using the data 
quality files and the WFIXUP task. For each galaxy, the align- 
ment of the images obtained with the same filter was checked 
by comparing the centroids of stars in the field of view. The im- 
ages were aligned to an accuracy of a few hundredths of a pixel 
using IMSHIFT and knowledge of the off'sets. They were then 
combined with IMCOMBINE. We verified that the alignment 

^ Imaging Reduction and Analysis Facilities (IRAF) is distributed by 
the National Optical Astronomy Observatories which are operated by 
the Association of Universities for Research in Astronomy (AURA) un- 
der cooperative agreement with the National Science Foundation. 



and combination did not introduce a significant blurring of the 
data. To this aim, we measured the FWHM of the Gaussian fit- 
ting to the field stars in the original and combined frames. We 
found that they did not change to within a few percent. In com- 
bining the images, pixels deviating by more than three times the 
local standard deviation - calculated from the combined eff'ect 
of Poisson and read-out noise - were flagged as cosmic rays and 
rejected. The residual cosmic rays and bad pixels were corrected 
by manually editing the resulting image with IMEDIT. 

The ACS images were calibrated using the CalACS reduc- 
tion pipeline. Reduction steps include bias subtraction, dark 
current subtraction, flat-fielding correction, and correction for 
geometric distortion with MULTIDRIZZLE task of IRAF as 
described in detail in ACS instrument and data handbooks 
(iPavlovskv et"aDl2004l2006h . The images obtained in the same 
filter were aligned by comparing the centroids of stars in the field 
of view and then combined, rejecting cosmic rays in the process. 
Residual cosmic ray events and hot pixels were removed using 
the LACOSJMA procedure (van Dokkum 2001). 

The sky level in WFPC2 images was determined from ap- 
parently empty regions in the Wide Field chips and subtracted 
from the PC frame after appropriate scaling. In ACS images, 
it was determined from regions free of sources at the edge 
of the field of view and then subtracted. The ACS/F330W 
and WFPC2/F336W passbands approximate Johnson-Cousins 
U band. The ACS/F475W, WFPC2/F555W and WFPC2/F8 14W 
passbands are similar to Johnson-Cousins B, V, and / bands, re- 
spectively. The flux cali bration to the Vega m agnitude system 
was perfor med following iHoltzrnan et al.l (1199 5) for the WFPC2 
images, and lSirianni et all (l2005l) for the ACS images. 

4. Unsharp images of nuclear disks 

To first gauge the structure and extent of the NSDs in all the 
available images of the sample galaxies, we constructed the 
unsharp-masked imag e of each f r ame u sing an identical proce- 
dure described by Piz zeUa et al.l (l2002h . The ACS/HRC image 
was resampled to match the same spatial scale as the WFPC2/PC 
images. Each image was divided by itself after convolution by 
a circular Gaussian of width cr = 6 pixels, corresponding to 
a'28 and a.' 29 in the WFPC2/PC and ACS/WFC images, respec- 
tively (Fig.[T]). This procedure enhanced any surface-brightness 



4 



L. Morelli et aL: Multiband photometric decomposition of nuclear stellar disks 




F336W F475W F555W " F814W 



Fig. 1. Unsharp masking of the HST archival images of NGC 4458 (ACS/F475W, WFPC2/F555W, and WFPC2/F814W; top panels), 
NGC 4478 (ACS/F330W, ACS/F475W, WFPC2/F555W, and WFPC2/F814W; middle panels), and NGC 4570 (WFPC2/F336W, 
ACS/F475W, WFPC2/F555W, and WFPC2/F814W; bottom panels). The orientation is specified by the arrow indicating north and 
the segment indicating east in the lower-left corner of the ACS/F475W image. The size of the plotted region is about 15'' x 15''. 



fluctuation and non-circular structure extending over a spatial 
region comparable to the cr of the smoothing Gaussian. Each 
galaxy nucleus clearly contains a highly elongated structure in 
all the images. The location, orientation, and size of each struc- 
ture is similar in all the observed passbands. These nuclear 
structures are associated with a central increase in ellipticity 
measured by performing an isophotal a nalysis using the IRAF 
task ELLIPSE (Figs. [2] - IH see also Scor za & van den Bosch 
[l998, Morelli et al. 2004). They are not artifacts of the unsharp- 
masking procedure as discussed by Morelli et al. (200 4|) for 
NGC 4458 and NGC 4478 and by .Scorza & van den BoschI 
(Il998h for NGC 4570. 

The ACS/F475W unsharped-masked image of NGC 4570 
exhibits two symmetric bright spots located along the major axis 
of the NSD at about 2" from the ce nter. They were identified by 
Ivan den Bosch & EmsellemI (Il998h with the signature of the nu- 
clear ring located at the inner Lindblad resonance in the equato- 
rial plane of the galaxy. This ring is seen edge-on. The two spots 
are barely visible in the WFPC2/F336W and WFPC2/F555W 
images and almost missing in the WFPC2/F 814W image. 
According to van de n Bosch & EmsellemI (Il998h . the blue col- 
ors of these features correspond to a rather young stellar pop- 
ulation with an age of less than 2 Gyr and a metallicity close 
to solar. The width of the filtering Gaussian we adopted in the 



unsharp-masking process is eff'ective in detecting the NSD, but 
does not enhance the structure of the ring located at the ultra- 
harmonic resonance nor the inner cutoff' to the main disk (for 
this purpose, see the unsharp-masked images obtained with dif- 
ferent filtering Gaussians that Ivan den Bosch & Emsellemll 19981 
displayed in their Fig. 1). 

In iMorelli et al.l (|2004|) . we derived isophotal profiles from 
the analysis of the WFPC2/F814W images of NGC 4458 and 
NGC 4478. Here, we present the isophotal profiles derived from 
the remaining images of the sample galaxies adopting the same 
method. We first masked foreground stars and then fitted el- 
lipses to the isophotes. We allowed the centers of the ellipses 
to vary, to test whether the galaxies were disturbed. Within the 
errors of the fits, we found no evidence of variations in the fit- 
ted center. The ellipse fitting was repeated with the ellipse cen- 
ters fixed. The signal-to-noise ratio of the ACS/F330W image of 
NGC 4478 did not allow us to measure a reliable radial profile 
of the a4 coefficient. This is necessary to derive the NSD param- 
eters, therefore we did not consider this image in the subsequent 
analysis. The resulting azimuthally averaged surface brightness, 
ellipticity, position angle, and the fourth cosine Fourier coeffi- 
cient (a4) profiles of NGC 4458, NGC 4478, and NGC 4570 
are presented in Figs. [2l [3l and IH respectively. In the inner- 
most about 1" of NGC 4458 and NGC 4570 and about 2" of 
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Fig. 2. Isophotal parameters of the nuclear region of NGC 4458 
as a function of the isophotal semi-major axis. We show the re- 
sults based on the analysis of the surface-brightness distribution 
measured in the ACS/F475W {left panels) and WFPC2/F555W 
{right panels) images, respectively, (a) Surface-brightness radial 
profiles of the galaxy after deconvolution {open black squares) 
and of the nuclear disk {solid blue line), (b) Radial profile of 
the fraction of total luminosity contributed by the nuclear disk. 
Radial profiles of the galaxy ellipticity (c), position angle (d), 
and fourth cosine Fourier coefficient (e) before {open black 
squares) and after (filled red circles) the subtraction of the best- 
fitting model for the nuclear disk. The results of th e photometric 
decomposition of the WFPC2/F814W image are in lMorelli et al.l 
(120041) . The vertical dash-dotted line indicates the radius inside 
which the deconvolution algorithm is not able to recover the ge- 
ometric parameters of the isophotes. 

NGC 4478, we measured positive values of the ^4 Fourier co- 
efficients, which describe the disky deviation of the isophotes 
from pure ellipses (lJedrzejewskilll987l) . These p hotometric fea- 
tures confirm the presence of a NSD (see also iKormendy et al.l 
I2009h . 

5. Structural parameters of nuclear disks 

After establishing the existence of NSDs in all passbands, 
we derived their photome tric parameters using the method 
of IScorza & Bended (Il995h as implemented by iMorelli et al.l 
(12004 . This method is based on the assumption that the isopho- 



Fig. 3. Same as in Fig. H but for NGC 4478. The results of the 
photometric decomp osition of the WFPC2/F814W image are in 
iMorelfi et al.1 (l2004l) . 



tal diskiness is the result of the superimposition of a spheroidal 
component (which is either a host elliptical galaxy or a bulge 
component) and an inclined exponential disk. The two compo- 
nents are assumed to have both perfectly elliptical isophotes with 
constant but diff'erent ellipticities. 

To properly derive the photometric parameters of the nu- 
clear disks, it is important to account for the HST point- spread 
function (PSF). The adopted PSF model was calculated with 
the TINYTIM package by taking into account the instrumen- 
tal setup and p osition of the nuclear disk on the given image 
dKrist&Hookl 119991) . To restore the images from the eff'ects 
of the P SF, they were deconvolved wi th the Richardson-Lucy 
method (lRichardsonlll972l: lLucvlll974l) by applying the IRAF 
task LUCY for a number of iterations between 3 and 6. A larger 
number of iterations was not found to improve the restoration of 
the ge ometrical param eters of the isophotes as discussed in de- 
tail by lMichardl 11996') with extensive experiments on the ground 
and HST-based images. The image noise was, instead, amplified 
and errors in the ellipticity and ^4 Fourier coefficients are en- 
larged. This eff'ect is illustrated in Fig.[5]for the WFPC2/F555W 
image of NGC 4458. The average surface brightness in the 
galaxy nucleus and the image noise are plotted as a function of 
the number of iterations performed with the deconvolution al- 
gorithm. The average surface brightness of the galaxy was mea- 
sured in circular apertures of radii 2, 3, and 6 pixels, respectively 




and centered on the nucleus. They correspond approximately to 
3,5, and 10 times the standard deviation of the best-fit Gaussian 
model of the PSF. The noise was estimated by measuring the 
standard deviation in the surface brightness level for a number of 
empty regions at the edges of the field of view. After 6 iterations, 
a negligible rise (< 2%) in the surface brightness is observed 
even in the innermost aperture, which is the most sensitive to 
the deconvolution eff'ect. In contrast, the noise level increases 
continuously at each iteration. This increases the uncertainties 
in the ^4 Fourier coefficients, which are critical for recovering 
the structural parameters o f the NSD with the met hod of pho- 
tometric decomposition by IScorza & Bended ([19951) . The radial 
profiles of the a4 Fourier coefliicients measured in the original 
WFPC2/F555W image of NGC 4458 as well as in the images re- 
sulting after 6 and 20 deconvolution iterations are shown in Fig. 
[51 The redistribution of the galaxy light caused by the deconvolu- 
tion is characterized by an enhancement of the galaxy diskiness 
in the inner V. The ^4 Fourier coefficients measured after 6 and 
20 iterations have the same radial trend and are consistent within 
the errors. However, the a4 errorbars increased by a factor of 3 
after 20 iterations. We found that after 40 Richardson-Lucy it- 
erations the scatter and errorbars of the ^4 Fourier coefficients 
are sufficiently large to blur the peak in the ^4 radial profile and 
significantly reduce the radial range along which the diskiness is 
observed. This prevents us obtaining accurate and reliable values 
of the NSD parameters. 



The rms tracking error is negligible with respect to the PSF 
FWHM (O'.'OT), thus it does not affect the image quality and the 
measurement of the structural parameters of the NSDs. To this 
aim, no correction for telescope jitter was necessary. 

ThelScorza & Bended (11995') method consists of the iterative 
subtraction of a thin disk model characterized by an exponential 
surface-brightness profile of central surface brightness yUo, scale- 
length h, axial ratio b/a, and position angle PA. The disk incli- 
nation was calculated to be / = arccos (b/a). The position angle 
of each nuclear disk was found to be constant within 1° for all 
the images of the same galaxy (Table [2l). Adopting this value, 
the remaining disk parameters are adjusted until the departures 
from perfect ellipses are minimized (i.e., ^24 0). 

iFerrarese et al ] (I2OO6I) measured in HST/ACS images the 
isophotal parameters of sample of early-type members of 
the Virgo cluster (including NGC 4458, NGC 4478, and 
NGC 4570). Inside r < 0.1 re, there is a trend for the mean ^4 to 
increase as galaxies become fainter. On large scales (OAr^ < r < 
Te), the isophotes are consistent with having a constant shape. In 
this radial range, most (but not all) of the observed galaxies ex- 
hibit almost elliptical isophotes. These results agree with the pre- 
vious findings based on HST photometry ([v an den Bo sch et al. 
"1994 lLauer etal.n i995- Verdoes Kleiin et al. il999i: (Rest et al. 
200 1| but see also Kormendy et al. 2009). We therefore decided 
to adjust the NSD parameters until the isophotes of the host 
spheroid had the same shape observed outside the inner region 
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Fig. 5. Deconvolution of the WFPC2/F555W image of NGC 4458 with the Richardson-Lucy algorithm. The left panels show (a) the 
average surface brightness per pixel in the galaxy nucleus, (b)its relative increase, and fcj the noise level of the image plotted as a 
function of the number of the iterations of the deconvolution process. The different symbols refer to the values measured in circular 
apertures with radii of 2 (green squares), 3 (red triangles), and 6 (blue circles) pixels centered on the galaxy nucleus, respectively. 
The relative increase in surface brightness at iteration 1 is calculated with respect to the original undeconvolved image. The noise 
level (filled circles) of the image was measured after each iteration in the same empty regions at the edges of the field of view. 
The right panels show the radial profiles of the fourth cosine Fourier coefficients measured in the (d) original image, and in the 
deconvolved image after (e) 6 and (f) 20 iterations. 



where ^4 > 0. We adopted ^4 ^ 0.01 for NGC 4458, ^4 ^ for 
NGC 4478, and ^4 ^ -0.01 for NGC 4570, respectively. 

For each disk model, the disk-free image of the galaxy was 
obtained by subtracting the disk model from the galaxy image. 
The isophotal analysis of the disk- free image was performed us- 
ing the IRAF task ELLIPSE. The photometric decomposition 
was performed independently for each image. The best-fit val- 
ues of yuo, h, and b/a and their 68% confidence-level err ors were 
derived in a similar way to those in lMorelli et al.l (|2004l) and are 
listed in Tab. [2j The comparison between the isophotal param- 
eters of the galaxies measured before and after the subtraction 
of the best-fit model of their nuclear disk are shown in Figs.O- 
m The surface-brightness radial profiles of the nuclear disks and 
the radial profile of the light fraction D/T that they contribute to 
the galaxy luminosity are also plotted. 

The a4 Fourier coefficients observed in the three galaxies 
drop to zero at small radii (< O'.'S), where the galaxy isophotes 
become rounder (e < 0.4) and perfectly elliptical (Figs. [2] - 
131). The NSD models have elliptical isophotes of higher ellip- 
ticity (e ^ 0.8; Table O. When they are subtracted from the 
galaxy images, the innermost residual isophotes becomes boxy 
(a4 < 0). This is an artefact of t he Richar d son-L ucy decon- 
volution method. As discussed by iMichardI (Il996l) . the algo- 
rithm fails to recover the geometrical properties of the galaxy 
cores. In particular, the restored ellipticity is underes timat ed 
and decreases to zer o towards the center (see also lLucyiri994) . 
Ivan den Bosch et aO (I1998i) place a lower limit of (YfS for re- 



covering relia ble geom etric parameters in HST/WFPC2 images, 
and Scorza & van den Bosch (1998) adopted this limit in their 
nuclear disk decomposition of NGC 4342 an d NGC 4570. We 
assume that the results obtained by van den Bosch et aD ([l998h 
are directly applicable to our case: because we study images ob- 
tained with similar or longer integration times, galaxies with 
shallower surface-brightness profiles and NSDs with equal or 
larger scale lengths. The decomposition results do not depend 
on the surface-brightness distribution in this region. Figure [6l 
shows the shape of the isophotes of the sample galaxies before 
and after the subtraction of the best-fit model of the NSD. The 
WFPC2/F555W images are considered here as an example. The 
diskiness of the galaxy isophotes disappears after the model sub- 
traction. The same behaviour is also observed in the other avail- 
able images. 

NGC 4458. To find the best-fit model of NGC 4458, we adjusted 
the disk parameters until the galaxy isophotes in the inner V 
have a4 ^ 0.01 instead of being perfectly elliptical (Fig.O. The 
scale length, inclination, and position angle of the nuclear disk 
measured in the ACS/F475W and WFPC2/F555W passbands are 
consistent within the errors. They a re also consist e nt wit hin the 
errors with the values measured by iMorelli et al.l (|2004|) in the 
WFPC2/F814W image (Fig. [7]). The mean values of the param- 
eters are (h) = 13 pc, (/) = 83°, and (PA) = 7°, respectively. It 
is the smallest and least luminous NSD (Lf814w = 3-10^ L©) 
of the sample. Its light contribution peaks at 0'.'3, where D/T 
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Table 2. Photometric parameters of the stellar nuclear disks. 



Galaxy 


Filter 




h 


i 


PA 


L 


Ref. 






[mag arcsec"^] 


["] [pc] 


n 


n 


[10" Lo] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 



NGC 4458 



NGC 4478 



NGC 4570 



F475W 
F555W 
F814W 
F475W 
F555W 
F814W 
F336W 
F475W 
F555W 
F555W 
F814W 




84^3 

78+}2 
78!l? 
80!|g 

78+12 
83!^ 



811 



7+1 
7+1 
7+1 
153+1 
153+1 
153+1 
158+1 
158+1 
158+1 

158+1 



2J+3.3 

3.3+^-^ 
3.0+^-^ 

17.0!^:^ 

4.5 

13.7__ 
1 9 q+i6h 
^^•^-5.0 
19 

99 4+22.0 



NOTES: Col.(l): Galaxy name. Col. (2): Filter name. Col. (3): Disk central surface-brightness in the Vega magnitude system. Col. (4): Disk scale 
length. A distance of 12.6 Mpc was adopted for all the sample galaxies. Col.(5): Disk inclinati on Col.(6) Disk position angl e. Col.(7): Disk 
luminosity. Col. (8): Reference paper for the photometric decomposition. 1 = this work, 2 - iMorelli et al I (120041) . 3 - iScorza & van den BoschI 
([T998h . 
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Fig. 6. Contour plots of the WFPC2/F555W deconvolved images 
of NGC 4458, NGC 4478, and NGC 4570 before (left panels) 
and after (right panels) the nuclear disk subtraction. Scales are 
in arcseconds, and orientations are the same as in Fig.[T] 



ranges from about 0.4 (ACS/F475W) to 0.2 (WFPC2/F555W) 
depending on the passband (Fig. [2]). 

NGC 4478. The scale length, inclination, and position angle of 
the nuclear disk in the ACS/F475W and WFPC2/F555W pass- 
bands are consistent wi t hin th e errors with the WFPC2/F814W 
values by iMorelli et al.l (|2004 Fig. (T]). The mean values of the 
parameters are (h) = 37 pc, </) = 78°, and (PA) = 153°, respec- 
tively. The nuclear disk contributes a maximum D/T = 0.2 in 
the radial range between 0'.'2 and (YfS from the center (Fig. [3]). It 
is one of the largest and more luminous NSDs (Lfshw = 2-10^ 
Lq) known to date. 

NGC 4570. To find the best-fit model of NGC 4570, we ad- 
justed the disk parameters until the galaxy isophotes in the in- 
ner V have ^4 ^ -0.01 instead of being perfectly elliptical 
(Fig. |4]). In this way, we account for the isophotal shape of 
the galaxy bulge, which is characterized by an increasing boxi- 
ness to a4 ^ -0.01 going from 3'' to V\ The central surface- 
brightness, scale length, and inclination of the NSD were al- 
ready obtained from the analysis of the WFPC2/F555W image 
by Ivan den BoschI (Il998l) . We also decided to perform the pho- 
tometric decomposition using our method too. Homogeneous 
decompositions and errors in the fitted parameters are required 
to properly compare the structural parameters of the NSD in 
the diff^er ent passbands. Our res ults agree within the errors with 
those by Ivan den Bosc 3 (119981) . Moreover, the averaged scale 
length, inclination, and position angle remain constant within the 
errors in all the observed passbands (Fig.|7l). Their mean values 
are (h) = 18 pc, </) = 81°, and (PA) = 158°, respectively. The 
disk luminosity ranges from 5 x 10^ (WFPC2/F336W) to 2 x 10'^ 
Lo (WFPC2/F814W). Its maximum D/T ^ 0.5 is observed at 
about a'3 (Fig.H. 

The trend of the central surface brightness in the diff'er- 
ent passbands is similar for all the measured NSDs. The value 
of the central surface brightness is constant within the er- 
rors in the ACS/F475W and WFPC2/F555W passbands and 
decreases by 1.6 mag arcsec"^ in the WFPC2/F814W band. 
The WFPC2/F336W image was available only for NGC 4570. 
The central value of the surface brightness in this filter is 
about 1.3 mag arcsec"^ greater than that in ACS/F475W pass- 
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Table 3. Colors of the NSDs and their host spheroids 



Color 


NGC 4458 


NGC 4478 


NGC 4570 




NSD Spheroid 


NSD Spheroid 


NSD Spheroid 


U -B 






1.1 0.8 


B-V 


0.6 1.0 


0.8 1.1 


0.9 1.0 


V-I 


1.4 1.1 


1.2 1.1 


1.5 1.2 



NGC 4458 
) _ NGC 4478 
_A NGC 4570 



NOTES: The conversion to the Johnson system has been calculated as 
U-B= [F336W-F475W]-0.21;5-y = [F475W - F555W] + 0.36; 
y - / = [F555W - F814W] - 0.05. 



band. In Table [3l we report the colors NSDs and of the 
their host spheroids. The colors of the spheroid were mea- 
sured along the disk major axis at a distance of 10 scale 
lengths on both sides of the center, where the contribution 
of the NSD to the total light is negligible (Figs. [2] - 131). 
The conversion to the Johnson system was calculated using 
SYNPHOT in IRAK Since this correction depends on the spec- 
tral energy distribution of the object, it was calculated using 
the iKinn ey et al. (1996) spectral templates. The resulting col- 
ors are similar to those measured in galaxi es hosting a NSD 
(NGC 4128, NGC 4621, and NGC 530 8, iKrainovic & Ja^ 
120041: NGC 4486A. lKormendv et al.ll200l ) The large uncertain- 
ties in the central surface brightnesses are intrinsic to the photo- 
metric decomposition method. This prevented us from perform- 
ing color analysis as a consistency check with the measuremen t 
of age and metallici ty performed by van den Bos ch et a (119981) . 
iMorelli et aP (l2004l) . and Krajnovic & Jaffe ( 2004|)and based on 
the line- strength indices. 

The scale length, inclination, and position angle derived for 
each NSD in the different filters are constant within the errors, 
independently of their values and the properties of their host 
galaxy. The nuclear disk has the same position angle as the inner 
part of the host spheroid. 

6. Discussion and conclusions 

We have investigated the photometric properties of the NSDs 
of three early-type galaxies in the Virgo cluster, NGC 4458, 
NGC 4478, and NGC 4570 by analyzing the WFPC2 and ACS 
images of their nuclei available in the HST Science Archive. 
The images were unsharp masked and their visual inspection 
suggested that the size, orientation, and location of each NSD 
was independent of the observed passband. To verify this, we 
derived the central surface brightness, scale length, inclination, 
and position angle of each NSD in all the pa ssbands by apply- 
ing the photometric decomposition method of IScorza & Bender 
ii995^. Some of the images had already been analyzed by 
& den Bosch et al. (1998, WFPC2/F555W for NGC 4570) 
and IMorelli et al.1 (120041 WFPC2/F814W for NGC 4458 and 
N GC 4478). Since we adopted an identical procedure to that 
of IMorelli etaP (l2004l) . we reUed on their WFPC2/F814W de- 
composition of NGC 4458 and NGC 4478 and repeated the 
WFPC2/F555W analysis of NGC 4570 to build homogeneous 
set of photometric parameters and corresponding errors for all 
the NSDs in all the passbands. 

The structural parameters (i.e., scale length, inclination, and 
position angle) of each NSD are constant within the errors in 
all the observed passbands, independently of their values and 
the properties of their host spheroid. This implies an absence of 
color gradients in nuclear disks and can be used to constrain their 



. 
90 - 



F475W F555W 
HST photometric band 



Fig. 7. Photometric parameters of the nuclear stellar disks of 
NGC 4458 (blue circles), NGC 4478 (green squares), and 
NGC 4570 (red triangles) in the WFPC2/F336W, ACS/F475W, 
WFPC2/F555W, and WFPC2/F814W passbands. 



star formation processes. NSDs are considered to be the end re- 
sult of the formation of stars from gas of external (i.e., accreted 
from galaxy neighborhoods or captured during a merger with a 
gas rich-companion; e.g., NGC 4458, NGC 4478 IMorelli et al.1 
2004) or internal origin (i.e., conveyed from oute r galactic re- 
gions by a bar; e.g., NGC 4570 van den Bosch & Emsellernl 
LI 998) piled up in the galaxy nucleus. Thus independent of the 
way in which the gas was originally accumulated towards the 
center, the dissipational formation process has to be invoked to 
account for the dynamically cold structure of the nuclear disks. 
This is hard to explain if NSDs assembled from already formed 
stars. We interpret the absence of color gradients in NDSs as the 
signature of a star formation event that occurred homogeneously 
over the entire disk. Although our samp le is adm ittedly very 
small, an inside-out formation scenario unoz-Mateos et al. 
2007) seems to be ruled out in NSDs given that this would pro- 
duce radial population gradients that are not observed. 

Acknowledgements. We would like to thank Enzo Brocato, Valentin Ivanov, 
Jairo Mendez-Abreu and Lodovico Coccato for their useful suggestions. LM 
is supported by grant CPDR06 1795/06 from Padua University. MC gratefully 
acknowledges the European Southern Observatory (ESO) for support via ESO 
Studentship program at the ESO Research Facilities in Santiago. This research 
has been made possible by support from Padua University through grants 
CPDA068415/06, CPDA089220/08, and CPDR095001/09 and from Istituto 
Nazionale di Astrofisica (INAF) through grant PRIN2005/32. This research 



10 



L. Morelli et al.: Multiband photometric decomposition of nuclear stellar disks 



has made use of the Lyon Extragalactic Database (LEDA) and NASA/IPAC 
Extragalactic Database (NED). 



References 

Athanassoula, E., Lambert, J. C, & Dehnen, W. 2005, MNRAS, 363, 496 
Baggett, S., McMaster, M., Biretta, J., et al. 2002, WFPC2 Data Handbook, 

Version 4.0, STScI, Baltimore 
Balcells, M., Graham, A. W., & Peletier, R. F. 2007, ApJ, 665, 1084 
Barnes, J. E. & Hernquist, L. 1996, ApJ, 471, 115 

Bertola, F. & Corsini, E. M. 1999, in Galaxy Interactions at Low and High 
Redshift, ed. J. E. Barnes & D. B. Sanders, lAU Symp. 186, ASP, San 
Francisco, 149 

Bertola, E, Corsini, E. M., Vega Beltran, J. C., et al. 1999, ApJ, 519, L127 
CaroUo, C. M., Danziger, I. J., Rich, R. M., & Chen, X. 1997, ApJ, 491, 545 
de Vaucouleurs, G., de Vaucouleurs, A., Corwin, H. G., et al. 1991, Third 

Reference Catalogue of Bright Galaxies, Berlin, Springer (RC3) 
Dotti, M., Colpi, M., Haardt, R, & Mayer, L. 2007, MNRAS, 379, 956 
Eliche-Moral, M. C, Gonzalez-Garcia, A. C, Balcells, M., et al. 2009, in 

Hunting for the Dark: The Hidden Side of Galaxy Formation, ed. V. P. 

Debattista & C. C. Popescu, AIP Conf Ser. in press (arXiv: 09 11.5624) 
Emsellem, E., Bacon, R., Monnet, G., & Poulain, R 1996, A&A, 312, 777 
Emsellem, E., Cappellari, M., Krajnovic, D., et al. 2007, MNRAS, 379, 401 
Emsellem, E., Monnet, G., Bacon, R., & Nieto, J.-L. 1994, A&A, 285, 739 
Falcon-Barroso, J., Bacon, R., Bureau, M., et al. 2006, MNRAS, 369, 529 
Ferrarese, L., Cote, P., Jordan, A., et al. 2006, ApJS, 164, 334 
Halliday, C, Davies, R. L., Kuntschner, H., et al. 2001, MNRAS, 326, 473 
Holtzman, J. A., Burrows, C. J., Casertano, S., et al. 1995, PASP, 107, 1065 
Hopkins, R R & Quataert, E. 2010, ApJ, submitted (arXiv: 1002. 1079) 
Jedrzejewski, R. I. 1987, MNRAS, 226, 747 
Kinney, A. L., Calzetti, D., Bohlin, R. C, et al. 1996, ApJ, 467, 38 
Kormendy, J. & Bender, R. 1996, ApJ, 464, LI 19 
Kormendy, J., Bender, R., Richstone, D., et al. 1996, ApJ, 459, L57 
Kormendy, J., Dressier, A., Byun, Y. I., et al. 1994, in Dwarf Galaxies, ed. 

G. Meylan & R Prugniel, ESQ Conf Ser. 49, ESQ, Garching, 147 
Kormendy, J., Fisher, D. B., CorneU, M. E., & Bender, R. 2009, ApJS, 182, 216 
Kormendy, J., Gebhardt, K., Fisher, D. B., et al. 2005, AJ, 129, 2636 
Krajnovic, D. & Jaffe, W. 2004, A&A, 428, 877 

Krist, J. & Hook, R. 1999, STIS Instrument Handbook Version 4.0, STScI, 
Baltimore 

Kuntschner, H., Emsellem, E., Bacon, R., et al. 2006, MNRAS, 369, 497 

Lauer, T. R., Ajhar, E. A., Byun, Y.-L, et al. 1995, AJ, 110, 2622 

Ledo, H. R., Sarzi, M., Dotti, M., Khochfar, S., & MoreUi, L. 2010, MNRAS, 

accepted (arXiv: 1005.2524) 
Lucy, L. B. 1974, AJ, 79, 745 

Lucy, L. B. 1994, in The Restoration of HST Images and Spectra - II, ed. 

R. J. Hanisch & R. L. White, STScI, Baltimore, 79 
Maraston, C. 1998, MNRAS, 300, 872 

McDermid, R. M., Emsellem, E., Shapiro, K. L., et al. 2006, MNRAS, 373, 906 
McMaster, M., Biretta, J., & et al. 2008, WFPC2 Instrument Handbook, Version 

10.0, STScI, Baltimore 
Michard, R. 1996, A&AS, 117, 583 

MoreUi, L., Halliday, C, Corsini, E. M., et al. 2004, MNRAS, 354, 753 
MoreUi, L., Pompei, E., PizzeUa, A., et al. 2008, MNRAS, 389, 341 
Munoz-Mateos, J. C, GU de Paz, A., Boissier, S., et al. 2007, ApJ, 658, 1006 
Pavlovsky, C, Biretta, J., Bohlin, R., et al. 2004, ACS Instrument Handbook, 

Version 5.0, STScI, Baltimore 
Pavlovsky, C, Koekemoer, A., Mack, J., et al. 2006, ACS Data Handbook, 

Version 5.0, STScI, Baltimore 
Peletier, R. R, Falcon-Barroso, J., Bacon, R., et al. 2007, MNRAS, 379, 445 
PizzeUa, A., Corsini, E. M., MoreUi, L., et al. 2002, ApJ, 573, 131 
Ravindranath, S., Ho, L. C, Peng, C. Y., Filippenko, A. V, & Sargent, W. L. W. 

2001, AJ, 122, 653 
Rest, A., van den Bosch, R C, Jaffe, W., et al. 2001, AJ, 121, 2431 
Richardson, W., H. 1972, 'Optical Society of American Journal, 62, 476' 
Rix, H. & White, S. D. M. 1990, ApJ, 362, 52 

Sandage, A. & Tammann, G. A. 1981, A Revised Shapley-Ames Catalog of 

Bright Galaxies, ed. Carnegie Institution, Washington D.C. 
Scorza, C. & Bender, R. 1995, A&A, 293, 20 
Scorza, C. & van den Bosch, R C. 1998, MNRAS, 300, 469 
Seifert, W. & Scorza, C. 1996, A&A, 310, 75 

Seth, A. C, Dalcanton, J. J., Hodge, R W, & Debattista, V R 2006, AJ, 132, 
2539 

Simoes Lopes, R. D., Storchi-Bergmann, T., de Fatima Saraiva, M., & Martini, 

R 2007, ApJ, 655,718 
Sirianni, M., Jee, M. J., Benitez, N., et al. 2005, PASP, 117, 1049 
Thomas, D., Maraston, C, & Bender, R. 2003, MNRAS, 339, 897 



Tran, H. D., Tsvetanov, Z., Ford, H. C, et al. 2001, AJ, 121, 2928 

TuUy, R. B. 1988, Science, 242, 310 

van den Bosch, F. C. 1998, ApJ, 507, 601 

van den Bosch, F. C. & Emsellem, E. 1998, MNRAS, 298, 267 

van den Bosch, F. C, Ferrarese, L., Jaffe, W, Ford, H. C, & O'Connell, R. W. 

1994, AJ, 108, 1579 
van den Bosch, R C, Jaffe, W., & van der Marel, R. R 1998, MNRAS, 293, 343 
van Dokkum, R G. 2001, PASP, 113, 1420 

Verdoes Kleijn, G. A., Baum, S. A., de Zeeuw, R T., & O'Dea, C. R 1999, AJ, 

118, 2592 
Worthey, G. 1994, ApJS, 95, 107 



